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The orbital motion of inspiralling and coalescing black hole binaries can be investigated using 
a variety of approximation schemes and numerical methods within general relativity: post- 
Newtonian expansions, black hole perturbation theory, numerical relativity, and the effective- 
one-body formalism. We review two recent comparisons of the predictions from these various 
techniques. Both comparisons rely on the calculation of a coordinate invariant relation, in 
the case of non-spinning binary black holes on quasi-circular orbits. All methods are shown 
to agree very well in their common domain of validity. 

1 Introduction 

The detection and analysis of the gravitational radiation from black hole binaries by the ground- 
based LIGO/ Virgo observatories, and future space-based antennas, requires very accurate theo- 
retical predictions for use as gravitational wave templates. The orbital motion of such compact 
binary systems can be analyzed using multiple approximation schemes and numerical meth- 
ods in general relativity: post-Newtonian (FN) expansions, black hole perturbation theory, the 
effective-one-body (EGB) formalism, and numerical relativity (NR). It is crucial to compare the 
predictions from these various techniques for several reasons: such comparisons (i) provide inde- 
pendent consistency checks of the validity of the various calculations, (ii) they help to delineate 
the respective domains of validity of each method, and (iii) they can inform the development of 
a universal semi-analytical model of the binary dynamics and gravitational wave emission. In 
this paper, we shall summarize the main results of two such recent comparisons, both related to 
the local orbital dynamics of non-spinning black hole binaries on quasi-circular orbits. 

2 Redshift Observable 

Gur first comparison is concerned with the relativistic motion of compact binary systems within 
black hole perturbation theory and the FN approximation. Consider two non-spinning black 
holes with masses mi and m-2, moving on an exactly circular orbit with angular frequency Q^p. 
The dissipative effects associated with the emission of gravitational radiation are neglected, 
which is formalized by assuming the existence of a helical Killing vector (HKV) field /c". The 
4-velocity of the "particle" mi is necessarily tangent to the HKV evaluated at that location; 
hence uf = U kf. The scalar [/ is a constant of the motion associated with the helical symmetry. 
It also measures the gravitational redshift of light rays emitted from mi, and received at large 
distance, along the helical symmetry axis perpendicular to the orbital plane;'^we shall henceforth 



refer to U as the "redshift observable". Being coordinate invariant, the relation U{^^p) provides 
a handy testbed to compare the predictions from the two approximation schemes. 

For an extreme mass ratio black hole binary, such that mi <C m2, the redshift observable 
U{yi^\mi,m2) is conveniently expanded in powers of the mass ratio q = mi/m2, according to 

U = Usch^ + qUGSF + 0{q^). (1) 

All coefficients in the expansion ([T]) are functions of the dimensionless coordinate invariant PN 
parameter y = (7712^^,^)^''^. The result for a test mass in circular orbit around a Schwarzschild 
black hole of mass m2 is known in closed form as C/schw = (1 ~ 3y)~^^^. The invariant relation 
f^GSF(y) encoding the first order mass ratio correction has been computed numerically, with 
high precision. EE This gravitational self-force (GSF) effect has also been computed analytically 
up to high PN orders. HH! The post-Newtonian expansion of f/csF is of the form 

^7GSF = ^afc/+i+lny^/3fe/+i + ... , (2) 

where the coefficients and /3k are pure numbers, and the dots stand for terms involving 
powers of logarithms {InyY, with p ^ 2, which are expected not to occur before the very high 
7PN order. ^2 The Newtonian, IPN, 2PN and 3PN polynomial coefficients {00,01,02,03} were 
determined analytically, EE as well as the leading-order 4PN and next-to-leading order 5PN 
logarithmic coefficients {P^, 13^} Their values are reported in the left panel of Table [l] 

Table 1: The analytically determined post-Newtonian coefficients Ok and fik (left panel), and the numerically 
determined values of higher-order PN coefficients, based on a fit to the GSF data (right panel). The uncertainty 
in the last digit is indicated in parenthesis. 
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Making use of the known results for the coefficients {aoi ^i, 02, /34, /^s}, a fit to the GSF 
data for UcsFi^ip) gave the numerical estimate a|* = —27.6879035(4) for the 3PN coefficient, SI 
to be compared with the exact value 03 = — 27.6879026 The results are in agreement 

with nine significant digits, at the 2a level. This provides a strong and independent test of the 
validity of both calculations, which rely on very different regularization schemes to subtract the 
divergent self-fields of point particles (mode-sum regularization in the self-force, and dimensional 
regularization in PN theory). By fitting the accurate GSF data to a PN model of the form 
([2]), now taking into account all known PN coefficients, including the exact value of the 3PN 
coefficient 03, the values of previously unknown PN coefficients Ok and /3k were measured, up to 
the very high 6PN order .1^1 These are reported in the right panel of Table [l] Notice in particular 
how the 4PN and 5PN coefficients 04 and as could be determined with high precision. 

Figure [1] shows the exact results for Uqsf{^ip), as computed within the self-force, as well as 
the successive truncated PN series up to 6PN order, based on the analytically and numerically 
determined PN coefficients summarized in Table [Tj This comparison illustrates the complemen- 
tarity of the two approximation schemes: previous knowledge of analytically determined "low" 
order PN coefficients allows to extract from the accurate GSF data information about higher 
order PN effects, which otherwise would likely remain inaccessible to standard PN calculations. 
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Fi gure 1: The gravitational self-force contribution f/csF to the redshift observable f/, as a function of rn = ra^jy, 
a coordinate invariant measure of the orbital separation. Notice that rn = 6m2 corresponds to the very relativistic 
innermost stable circular orbit (ISCO) of a test-mass orbiting a Schwarzschild black hole of mass mi. 

3 Periastron Advance 

As long as the radiation-reaction time scale is much longer than the typical orbital time scale, the 
motion of two non-spinning black holes on a generic eccentric orbit depends on two independent 
frequencies: the radial frequency (or mean motion) il^ = 27r/P, where P is the radial period, 
i.e. the time interval between two successive periastron passages, and the periastron precession 
frequency A$/P, where A$/(27r) = K — \\^ the fractional advance of the periastron per radial 
period. In the zero eccentricity limit, the relation between the circular orbit frequency and 
K = n^p/V^r is coordinate invariant; it can thus be used as a convenient reference for comparison. 

The invariant relation K{rt^) has been computed at the 3PN accuracy in FN theory, El at 
first order in perturbation theory,!^ and in the FOB formalism.^ This genuine general relativistic 
effect has also recently been measured for the first time in fully non-linear NR simulations.'^ 
Le Tiec et aZ.!^ considerably improved upon the accuracy of this initial measurement. Making 
use of new and longer simulations of the late stage of the inspiral of non-spinning black hole 
binaries with mass ratios q = 1,2/3, 1/3, 1/5, 1/6, and 1/8,^ they measured K with a relative 
uncertainty ~ 0.1 — 1%. This accuracy made possible an extensive comparison which, for the 
first time, (i) encompassed all the analytical and numerical methods currently available, and (ii) 
focused on the orbital dynamics of the binary, rather than the asymptotic waveform. 

Figure [2] shows the invariant relation K{Q.ip) for binary black holes with mass ratios (? = 1 
(left panel) , and g = 1/8 (right panel) , as computed in NR (in cyan) , FN theory (red) , and the 
FOB formalism (yellow). For comparable masses (e.g. g = 1 or 2/3), the 3PN prediction is in 
good agreement with the exact result from NR (to better than 1%). However, as expected, it 
performs less well when q — )■ O.'^ The FOB (3PN) prediction, on the other hand, is in very good 
agreement with the NR data over the entire range of frequencies and mass-ratios considered. 

Also shown in Fig.[2]are the predictions for a test mass in circular orbit around a Schwarzschild 
black hole (green), and the inclusion of the GSF (magenta and blue). While perturbative self- 
force calculations are commonly formulated as expansions in powers of the usual mass ratio q 
[see e.g. Fq. Q], FN expansions naturally involve the symmetric mass ratio v = mim2/m^, 
where m = nii + m2 is the total mass of the binary. Since at first order q = u + 0{u^), the GSF 
result for the periastron advance may as well be written in the "resummed" form 




(3) 




Figure 2: The periastron advance K — 1 + A$/(27r), as a function of the circular orbit frequency f2^, for black 
hole binaries with mass ratios 1 : 1 (left panel) and 1 : 8 (right panel). Notice that in the later case, A<E> reaches 
half an orbit per radial period for mO,^ ~ 0.03, corresponding to an orbital separation r ~ 10m. 



where all coefficients are functions of the dimensionless invariant PN parameter x = (mO,^)^/'^. 
The GSF correction Kqsf to the test-particle result i^schw = (1 — 6x) ' has recently been 
computed numerically.'^ Although the GSFg prediction [obtained by replacing u ^ q in Eq. ([s])] 
agrees with the exact result within a relative difference of magnitude ~ g^, as expected, the GSFj^ 
prediction ([s]) agrees remarkably well with the NR data for all mass ratios. This surprising result 
suggests that GSF calculations may very well find application in a broader range of physical 
problems than originally envisaged, including the modelling of intermediate mass ratio inspirals, 
a plausible source of gravitational waves for Advanced LIGO/Virgo. 
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